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We present an approach to electronic polarization in molecular solids treated as a set of quan- 
tum systems interacting classically. Individual molecules are dealt with rigorously as quantum- 
mechanical systems subject to classical external fields created by all other molecules and, possibly, 
external sources. Self-consistent equations are derived for induced dipoles and for atomic charges 
whose redistribution in external fields is given explicitly by an atom-atom polarizability tensor. 
Electronic polarization is studied in two representative organic molecular crystals, anthracene and 
perylenetetracarboxylic acid dianhydride (PTCDA), and contrasted to previous results for systems 
of polarizable points. The stabilization energies of the neutral lattice, of isolated anions and cations, 
and of cation-anion pairs are found. Charge redistribution on ions is included. The dielectric ten- 
sors of crystals are successfully related to gas-phase properties and provide consistency checks on 
polarization energies. The procedure is generally applicable to organic crystals in the limit of no 
intermolecular overlap. 



I. INTRODUCTION 

Rftcent advances in the preparation j-of ordered thin 
filmsEl and of organic molecular crystalsoa with reduced 
impurity levels have revived interest in organic electron- 
ics and the properties of charge carriers in these mate- 
rials. Organic molecular crystals are typically insulators 
with low dielectric constant, k ^ 3, and charges local- 
ized on molecules. Electronic polarization, an effect that 
is usually not considered in conventional inorganic semi- 
conductors, has a central role in the electronic propertieia 
of organic crystals, as-jdiscussed by Gutmann and Lyons,u 
Pope and Swenberg,El and Silinsh and Capek.El When a 
charge carrier is brought into a molecular solid, its field 
polarizes the surrounding molecules. Secondary polariza- 
tion fields created by polarized molecules contribute to 
the total self-consistent polarization cloud that surrounds 
each charged quasiparticle. Such a cloud is sometimes 
referred to as an electronic polaron, to distinguish from 
lattice relaxation effects. 

The overall energy relaxation of a positive and nega- 
tive charge carrier, P+ + P_, is an inmortant property of 
an organic material. Typical valuesoQ are in the range 
of 2 — 3 eV. The transport gap, i.e. the energy necessary 
to create a well-separated electron-hole pair is 

Et = I-A-P+-P^ (1) 

where / and A are the gas-phase ionization potential and 
electron affinity of the molecule. P+ and P_ contain po- 
laronic contributions due to both intramolecular and lat- 
tice phonons that arc cstimatedQ to be ^ 10%. We focus 
on the large electronic components of P-|_ and P_ and 
from now on exclude polaronic effects. 

The appearance in Eq. of gas-phase properties is 
made possible by weak intermolecular forces and Van der 
Waals separations in organic molecular crystals. The 
solid-state environment is taken as a perturbation in 
molecular exciton theory. In contrast to inorganic semi- 
conductors, organic crystals can normally be approxi- 



mated as molecules with negligible overlap, and vanish- 
ing intermolecular overlap is the crucial approximation 
to our general development of electronic polarization en- 
ergies. The energy necessary to create an electron-hole 
pair at finite separation r defines the interaction potential 
V{r) < 0: 

Ppah.(r) =Pt + nr), (2) 

Nearest-neighbor ion pairs of large molecules obviously 
deviate from point charges. 

Mott and LittletonD first estimated the electronic po- 
larization energy in an atomic lattice by considering each 
atom as a polarizable point. The self-consistent treat- 
ment of point charges in lattices ofnpolarizaiblc points was 
subsequently developed by MunnJS Silinsh,l3 and cowork- 
ers. While satisfactory for atomic lattices, polarizable 
points for molecules -Completely neglect their structure. 
Polarization studies,Q'u primarily on the acene family, 
addressed molecular size to some extent by introduc- 
ing submolecules, whose choice is arbitrary, that carry 
a fraction of the total molecular polarizability. A recent 
applicationo to anthracene has several choices, the most 
elaborate. one being submolecules at all carbons, while 
a studyEj of perylenetetracarboxylic acid dianhydride 
(PTCDA) has 11 submolecules centered at rings and CO 
bonds. In effect, the quantum nature of molecules is ap- 
proximated by the microelectrostatics of submolecules. 

In large 7r-conjugated molecules subject to an exter- 
nal field, charge redistributes over distances comparable 
to the size of the molecule and generates large nonlinear 
optical responses. Such a flow of charge in molecules cre- 
ates secondary polarization fields that do not necessarily 
reduce to the field of a set of induced dipoles. Thus, rig- 
orous treatment of the polarization field requires analysis 
of the molecular charge distribution p(r). 

Another issue is electrostatic interactions present al- 
ready in the ground state of a condensed phase com- 
posed of neutral molecules, especially when they contain 
heteroatoms. These charges induce mutual polarization 
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in the surrounding molecules and contribute to the over- 
all stabilization energy of the solid. We estimate that 
the polarization contribution to the ground-state energy 
can be quite significant, reaching hundreds of meV per 
molecule. A notable exception is the acene family, whose 
TT-systems have approximate electron-hole symmetry and 
negligible partial charges, making the polarization energy 
small. 

In this paper we present an approach to electronic po- 
larization in molecular solids that allows for quantitative 
description of intramolecular charge redistribution. The 
crucial approximation is the neglect of intermolecular 
overlap. Zero overlap implies that coordinate space can 
be subdivided into non-overlapping regions, e.g. Wigner- 
Seitz cells, associated with individual molecules. Each 
molecule is then a quantum-mechanical system subject 
to external fields created by the crystal and, possibly, 
other external sources. The external fields are rigorously 
classical, so that quantum mechanics is needed at the 
intramolecular level only. 

Section II illustrates the idea and presents self- 
consistent equations. Section III describes an approx- 
imate discrete form of the equations that, we believe, 
strikes an optimum balance between accuracy and sim- 
plicity for practical use. In Sec. IV the equations are 
applied to the translationally-invariant lattice to find 
the polarization contribution to the binding energy. In 
Sec. V equations for the polarization energy of a system 
of ions embedded in a lattice are derived. The polariza- 
tion energy of isolated charge carriers in anthracene and 
PTCDA crystals is calculated, as well as the energy of 
various ion pairs. The dielectric tensors of anthracene 
and PTCDA are computed in Sec. VII, and their consis- 
tency with polarization energy calculations is established. 
We compare P_, and V{r) to submolecular results 
and comment briefly in Sec. VIII on applications to or- 
ganic solids with mobile charge carriers. 



II. SELF-CONSISTENT CHARGE 
DISTRIBUTION 



In the zero-overlap approximation, the self-consistent 
solution of the Schrodinger equation for a solid reduces to 
a product of wave functions of individual molecules. The 
minimum energy relative to gas-phase (nonintcracting) 
molecules or ions can be done in two steps. First, the 
change E{(j)) in the ground state energy of each molecule 
is found as a functional of the external electrostatic po- 
tential (/'(r). For simplicity we neglect any magnetic in- 
teractions. The ground state charge distribution p(r; (j>) 
also depends on <j) and determines the secondary polar- 
ization field created by the molecule. 

The total electrostatic potential at a point r within a 
molecule a is created by all other molecules b ^ a and, 
possibly, an applied field: 



(3) 



The prime at the sum excludes the term with b = a. The 
total energy of the solid is then 



a a<b'' 



(4) 



The second term compensates for double-counting the in- 
termolecular interactions in the first term. Equivalently, 



Etot - 



Cppi)p" 



(5) 



Minimization of i?tot with respect to (r) yields the self- 
consistent ground state energy of the molecular solid in 
the approximation of zero overlap. 

In principle, variation of Eq. (^) with respect to 6°- 
gives an equation for p°'{r), which together with Eq. (J^) 
forms a complete self-consistent system. A more efficient 
way, perhaps, is to use an appropriate quantum chemical 
procedure to find p°(r; 0") for each molecule a and iter- 
ate several times by updating ^''(r) using Eq. M). The 



]) and (|) for 
electrostatics. 



self-consistent problem defined by Eqs. ( 
charges and potentials is typical of classica 
The quantum part is limited to the charge distribution 

The index a can be dropped or restricted to a single 
unit cell when the translationally-invariant state of the 
lattice is of interest. This is the case, for example, in the 
calculation of the dielectric tensor .lid Otherwise, a finite 
number of molecules must be considered. A practical 
implementation of the procedure requires some form of 
discretization of the continuous functions cj)'^ (r) and p° (r) 
defined within the molecular volume. Certain trade-off 
between the accuracy and simplicity suitable for repet- 
itive quantum-chemical calculation is unavoidable. In 
the following section we develop a simple scheme, which 
captures intramolecular charge redistribution to a quan- 
titative accuracy. The procedure has been successfully 
implemented by us to calculate indices of refraction of 
anthracene and PTCDA.B 



III. MOLECULE IN NONUNIFORM FIELD 

In this section we omit the index a and consider a sin- 
gle molecule subject to an external potential <j){r). We 
note that charge redistribution gives a major contribu- 
tion to the polarizability of large conjugated molecules. 
This "major part" is not defined quantitatively, as there 
is no unique definition of atomic charges. The scheme we 
develop below separates molecular polarizability into two 
parts, the sum of which matches the actual molecular po- 
larizability a with the best value known from experiment 
or theory. 
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We use a semiempirical Hamiltonian because it pro- 
vides a natural way to represent an arbitrary external 
potential acting on a molecule. We define (j)i = (l){ri), 
the potential at the position of each atom r^. A site en- 
ergy is added to the diagonal matrix elements for the 
orthogonalized valence orhitals of atom i. We employ 
the INDO/S Hamiltonian jl3 which is known to approxi- 
mate molecular properties at only a tiny fraction of the 
cost of ab-initio calculations. Throughout the paper we 
use Lowdin charges pf, where i labels atoms in molecule 
a. The charges are defined as the sums of occupation 
numbers of orthogonalized orbitals of atom i. 

The corresponding contribution a'^ to the actual po- 
larizability a is clearly restricted to the molecular plane 
in conjugated molecules. We associate the difference 
a — a'~' between the actual and INDO/S polarizabilities 
with "atomic" contributions caused by the distortion of 
atomic orbitals in the field. Atomic contributions are 
small corrections to large in-plane polarizabilities. We 
note that any choice of Pi{F) leads to in-plane a'^ as a 
consequence of a discrete charge distribution. 

Atomic contributions to a can be described, as in 
atomic lattices, in terms of induced dipoles situated at 
the positions of atoms. Based on such an idea we pro- 
pose the following minimal scheme that is designed to 
capture both charge-redistribution and "atomic" parts 
of the molecular response to external fields. We describe 
the state of an iV-atom molecule by a set of 2N vari- 
ables. Pi and fii, which represent the partial charges and 
induced dipole moments of atoms. The same number of 
variables, <j>i — (/'(rj) and Fi ~ — V0(ri), describes the 
external field acting on a molecule. 

We denote by qi the deviation of partial charges from 
the ground-state values pf''^ of an isolated molecule or 
molecular ion, 



(0) 



(6) 



At small fields the energy of the molecule is quadratic in 
the distortion from equilibrium: 

ij i 



Here the positive-definite charge stiffness matrix H^^ de- 
scribes the increase in the internal energy of the molecule 
when the charge distribution deviates form its zero-field 
equilibrium; the tensor plays the same role for 

atomic dipoles. At given configuration {(pi, Fi} of the ex- 
ternal field, the minimum of the energy functional Eq. (j^) 
is achieved at 



p.=,f)-5:H, 



Hi = UiFi 



(8a) 
(8b) 



We see that cti is the polarizability for atom i. It 
does not necessarily reduce to a scalar, since atoms in 
a molecule have no rotational symmetry. We could also 
assume nonzero atomic dipoles iJ-f "' in the ground state of 
nonpolar molecules and so obtain symmetric equations, 
.(0) 



but we set fi] = in this paper. 

The energy of the molecule at the minimum is 



(9) 



Equivalently, 

E{(P, F) = Y, pfU^ - i E '^^n^,'^, -\Y. F^aF,. 

i ij i 

(10) 

The last two terms describe the energy relaxation of the 
molecule in the external field. The positive-definite sym- 
metric matrix H^ is the susceptibility with respect to site 
potentials (pi [cf. Eqs. (pa|) and (|lC 



H,; 



dpi 



•^3 / Q 



f d^E 

\d(j)^d(j)j 



(11) 



Hy deter- 



Partial derivatives are evaluated at (pi = 0. 
mines the charge redistribution among atoms in the ex- 
ternal potential. It is a natural extension of the similar 
quantity tt^ used in 7r-electron theoryJl3 and called the 
atom-atom polarizability. In our case, the total charge of 
all valence electrons is considered. Note that our defini- 
tion differs by a factor —1/2. 

Atom-atom polarizabilities H^ obey the condition 



= 0, 



(12) 



since the charge distribution in Eq. ( |8a| ) is invariant to 
an additive constant in all site potentials (pi. The zero- 
overlap approximation conserves charge at each molecule. 
Expansion to second order in (pi is sufficient for (pi < 1 eV. 
Eqs. (^) eliminate the need to solve repetitively the quan- 
tum problem for the molecule. We calculate N{N + 1) /2 
atom-atom polarizabilities Hy only once using Eq. (O) 
for the neutral molecule and for the cation and anion. 
Stronger perturbations may require re-evaluation of Hy 
at some intermediate (pi. 

The total induced moment of a molecule is 



E ^'^^P' 



fj-i 



(13) 



The molecular polarizability consists, therefore, of two 



terms, a — 4- a: 



-a/3 



(14) 



3 



where the Greek indices take the values x, y, and z. 

Equation (|lj) illustrates the advantages and limita- 
tions of partial atomic charges. With the aid of Ily, 
they provide a rigorous description of charge redistribu- 
tion. The assumption of polarizable points is atomic lat- 
tices is kept, however, through a — a — a'^ . We have 
corrections to INDO/S charges and distribute a propor- 
tionally to the numbers of valence electrons ni associated 
with iiidijddual atoms: Si = arii/ 'Y^rii. As in previous 
theory,El LJ a is an independent gas-phase input to the 
calculation. 

Table I Principal components of the molecular polarizabil- 
ities of anthracene and PTCDA; the long (L), medium (M), 
and normal (N) axes are fixed by D2h symmetry. 



pa 



Method 



aNN (A^) aMM (A^) ULL (A^) 



Anthrajcene 

Experimcnttil 15.2 25.6 35.2 

15.9 24.5 35.9 

B3LYP/6-311-F+G** 12.03 24.27 42.56 

INDO/S {a^) 24.05 41.52 

PTCDA 

B3LYP/6-31H-+G** 18.06 50.27 88.18 
INDO/S {a^) 50.84 84.54 

Table I summarizes results for anthracene and 
PTCDA. Density functional (B3LYP) reS|Blts have been 
obtained using the Gaussian 98 program.lla Theory and 
experiment are in reasonable agreement jJor. anthracene 
molecules when large basis sets, are usedJl^l-EI (such as 6- 
311-l-hG**). Dielectrifi-^atakl for crystalline PTCDA 
films are also consistently with calculated molecular po - 
larizabilities. The INDO/S results for a'^ from Eq. (|l|) 
are confined to in-plane components that represent 
charge redistribution according to liij in Eq. (|l4|). We 
also need a = a — a'-^. Unless otherwise indicated, we will 
use the B3LYP polarizabilities in Table I. Since they ex- 
ceed oP ^ atomic contributions increase the polarization 
compared to the "charges-only" choice of S = 0. We note 
that simple Hiickel theory often overestimates responses 
to applied fields and hence the amount of charge redis- 
tribution; in that case S may be negative and the atomic 
part reduces the polarization. Equations (H) hold for any 
a. 



IV. SELF-CONSISTENT EQUATIONS 

In the condensed phase the potential and field at the 
position of atom i of molecule a created by all other 
molecules h ^ a are 



(15a) 



b j 



bp 

3 ' 



(15b) 



where v{r) = 1/r, Va{r) = —dv/dr"", Vapiv) = 
d^v / dr"'dr^ . Summation over repeated Greek indices 
a, (3 = x,y,z is assumed. The vector rfj — — 
points to the atom of interest from atom j of molecule b. 
Here we have assumed no external sources for simplicity. 
Equations (^5|) together with Eqs. (^) form a complete 
self-consistent linear system for fif, (j>f, and Ff. 
The total polarization energy of the solid is [compare 

to Eq. m 



Et. 



E 



(16) 



After some algebra and using Eq. (||), this reduces to 



The derivation of this formula is simplified by replacing 
the dipoles by pairs of charges separated by small dis- 
tances, and taking the limit in the final expression. 

The total energy is a bilinear form of unperturbed 
charges p"^"'' and self-consistent potentials cf)'^. We can 
write it also in terms of self-consistent charges and 



dipoles ^f, defining the unperturbed potentials 

fields in Eqs. (15) by setting pf 

Using the identity 



,a(0) 



and 



) and (if = 0. 



Ep:^% 



EE-(r")A 

a^b ij 



aiO)^b 



Ep?^: 



the total energy becomes 

^^tot = iEE(^^^ 



q(0) 



(18) 



(19) 



We use this form of Etot in Sec. VI to treat ions in infinite 
lattices. r-inn 

Equation (p^ ) reduces to the previous resulttZfQD when 
all molecules are shrunk to points, with J2i P? = fo^' 
molecules and ±1 for ions. Since charge redistribution is 
no longer possible, we may set Ily = 0, = 0, and po- 
larizability 5; = a at positions of neutral molecules. The 
first correction to for finite molecules is an induced 
dipole at the center, which gives the same approxima- 
tion for a when combined with a. The potential (/'°^°'' in 
Eq. (|l|) is due to ions, and the first sum, which is now 
restricted to charged sites, becomes the source term Wq 
of Ref. |l^. The second term, over molecules, describes 
induced dipoles in the field of the ions and is the Wi 
term of Ref. The polarizability of ions is generally 
different from molecules, but is not required for finding 
P± in centrosymmetric lattices of point molecules, since 
the ion is at an inversion center. 
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The expression ( [T7| ) or (|T^) for the lattice polariza- 
tion energy is not restricted to equivalent molecules. In 
principle, each molecule a may have its own p. 



a(0) 



n? 



and af. In practice, there are several molecules per unit 
cell in organic molecular crystals. The translationally- 
invariant lattice of neutral molecules, the neutral lattice 
of the following section, reduces to atomic charges and 
potentials within a unit cell. Molecular ions in specified 
unit cells break translational symmetry and, as discussed 
in Sees. VI and VII, require different methods for finding 
-Etot- In the zero-overlap approximation, charge carriers 
are molecular ions in place of neutral molecules. The po- 
larization energy P± of a carrier is the energy difference 
between two extensive quantities, the lattice with the ion 
and the neutral lattice. 



therefore, reduces to a system of ANNc linear equations, 
where N(. is the number of molecules in a unit cell, and 
N is the number of atoms in a molecule. 



Madelung-type infinite sums-in Eqs. (15) can be eval- 
uated using Ewald's method.tS Special care has to be 
taken to treat complex lattices with many partial charges 
in a unit cell. For this purpose we introduce a fictitious 
uniform neutralizing background for each partial charge. 
Since the unit cell is neutral, these backgrounds cancel 
exactly. We define an auxiliary potential function V(r), 



V(r) = lim 

R^OQ 



E 



r-. 



2TTr' 

3Vr 



-( 



1/3 



(20) 



V. NEUTRAL LATTICE 

In this section we evaluate the polarization energy of 
the neutral lattice. The analogous quantity vanishes 
identically in the polarizable-point approach, since there 
are no fields or induced dipoles in the lattice until charges 
are introduced. The self-consistent Eqs. (||), ( p^ ) can be 
restricted to a single unit cell of volume Vc- The problem. 



V(r) 



3Vr 



erf(Gr) 



^'l-crf(G'|r-£|; 



e 



\r-i\ 



where the summation is over Mn lattice vectors £ falling 
within a sphere of radius R, with the term £ — missing. 
Subtracted is the potential of the uniform neutralizing 
spherical charge —Mn centered at the origin of coordi- 
nates. Centering of the neutralizing backgrounds at a 
common point in space is necessary for proper cancela- 
tion. 

Ewald's method gives 



TT cxp(-g74G2) 



E 



g2/4G2 



cos(gr), 



(21) 



where erf(a;) ~ (2/y^) dy exp{—y^) is the error func- 
tion. The second sum is over all reciprocal vectors g ^ 0, 
exp(jg£) = 1. Ewald's parameter G is arbitrary (the re- 
sult does not depend on its value); a reasonable choice 
is G = {tt'^ /vcY^^ . The function V(r) is regular within 
the central lattice cell £ = 0, including the point r = 0, 
and it is not periodic because of the missing term. The 
function V(r) + l/r is periodic in r. 

Equations ([l5| ) can be written in terms of V(r) and 
its derivatives VQ(r) = — 9V(r)/9r" and VQ^(r) = 
d'^V{r)/dr°'drf^: 

'^"-E'E^'(^^v?+^/3(r?f)/^f 

b 3 

+ E E ^("^i^'j + ^piOf^ (22a) 

b j 
b 3 

+ E E ^o.{vf)p) + V.Mr^')Mf • (22b) 

b j 

The sums over h are restricted to the central unit cell. In 
the primed sums the term 6 = a is excluded. The terms 
with V and its derivatives give contributions by charges 
and dipoles beyond the central cell. 



Equations ( p^ ) express 0° and F°; in terms of and 
/xj. Together with Eqs. (||) they form a complete linear 
system of ANNc equations, half of them vector. For ex- 
ample, NNc = 48 and 76, respectively, for anthracene 
and PTCDA lattices, which results in 192 and 304 scalar 
linear equations for these materials. 

The solution for the neutral lattice is further denoted 
as p°, /J", 4>i 1 and . These quantities are summed over 
molecules in Eq. ( |l7| ) or ([l9| ) to yield the (extensive) self- 
consistent energy of the neutral lattice. 



A. Anthracene and PTCDA 

Using the procedure described above we calculated 
polarization energy of anthracene and PTCDA crys- 
tals, which represent two major families of organic semi- 
conductors. Both materials are monoclinic with two 
molecules per unit cell. Both molecules have centers of 
inversion and are nearly planar in the crystal. PTCDA 
molecules are co-planar, up to a small tilt, and form lay- 
ers and stacks. In anthracene the angle between molecu- 
lar planes is significant. 

We used the X-pjy crystal structures for PTCDAEHI 
and for anthracenelij. The positions of hydrogens, not 
given accurately by X-ray, were AMl-optimized using 
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Gaussian Jli The gas-phase polarizabiUty, needed to de- 
termine the atomic correction a, is given in Table I. 

We obtain polarization energy of 330 meV per PTCDA 
molecule. This is two orders of magnitude greater than 
2.8 meV that we get for anthracene. The large polar- 
ization energy of the PTCDA lattice is caused by signifi- 
cant partial atomic charges in neutral PTCDA molecules, 
which are negligible in anthracene due to the approxi- 
mate electron-hole symmetry, as mentioned above. 




10 12 14 16 18 20 
Atom number 

Fig. 1 Partial charges in PTCDA molecule in crystal lat- 
tice (pi, solid line) and in gas phase {pf'\ dashed line) 

Figure 1 compares partial atomic charges of a PTCDA 
molecule in the gas phase and crystal. The inset ex- 
plains the atom numbering scheme. Only one half of 
the molecule is shown because of C^ symmetry. Charge 
redistribution yields excess positive charge on three hy- 
drogens whose partial charges roughly double. It is worth 
noticing that these hydrogens reside in positions that 
suggest the formation of incipient hydrogen bonds; the 
approximation of zero overlap excludes any covalent con- 
tribution. The distances from CH carbons to the nearest 
oxygen atoms in neighboring molecules are 3.338, 3.269, 
and 3.768 A for C5 — Oi', C14 — O19', and Cg — O19', re- 
spectively, while the corresponding C — H — O angles are 
143.8, 149.3, and 158.4° (see Fig. 2). 




Fig. 2 Arrangement of PTCDA molecules in a layer [pro- 
jection onto (102) crystalline plane]. Incipient hydrogen 



bonds (dotted lines) involve hydrogen atoms 5(a), 14(b), and 
6(c), according to the numbering. Crystalline coordinates 
here and in Table II conform to the notation of Ref. 20l 



VI. POLARIZATION ENERGY OF CHARGE 
CARRIERS 

We now consider a lattice with one or more neutral 
molecules replaced with molecular ions. To evaluate the 
energy, we solve self-consistent equations with the ions 
and subtract the polarization energy of the neutral lat- 
;h ion is described by n*°", 3*°", and p*""*-"-* 
^ = ±1), which may differ from the similar 



tice. 



quantities in the neutral molecule. While II-^" and 
are determined by semiempirical calculation, the atomic 
correction S*"" depends on a separate calculation for the 
ion. In this work, we use the same a for molecules and 
ions. 

It is useful to rewrite the self-consistent equations (H) 
and (|l^) in terms of the deviation from the neutral-lattice 
solution: 



Pi 



Pi. 



The equations for 6pf, 
and SF°' then read 



S<j,f = 0? - , (23a) 
5Ff^Ff~Fl (23b) 

6°, and each component of 



sp'} ^ f^* ^Y.^f 



~a3 r T-ia/? 



(24a) 
(24b) 



and 



-EE ^(r^')^/'? + Mrt'Wf, (25a) 

b j 

^-{r^'M + vo.0{vlf)5iif. (25b) 

b j 

The source term is zero except for ions, 

(fr* = ^pr^""^ - E An-"^°" ^ 5(fr(^). (26) 



zon(O) 



2on{0) 



.(0) 



An'°" 



TTion _ TT . . 



where Ap^ ' ' — Pj — p. 

and (pj is the potential at the ion's position in the lat- 
tice. In a hypothetical situation when a neutral molecule 
is replaced with a foreign molecule which, subject to 
the potential (j)^, has charge distribution p^, the source 
term is zero, and such a molecule will not disturb the 
translationally-invariant self-consistent solution. 

With ions present, the problem has no translational 
symmetry, and the number of self-consistent equations is 
infinite. We consider an imaginary cluster within an infi- 
nite lattice that includes all the unit cells within a certain 
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distance R from the origin. Some dozens of molecules 
are required for the cluster to resemble a sphere. We 
set Spf = = for the molecules outside the cluster, 
and solve the self-consistent Eqs. (^) and (|2^) for Spf, 
Sfif, S(l)f, and SFf within. This corresponds to an infi- 
nite lattice in which only the charges within the cluster 
are allowed to relax. Molecules outside retain the charge 
distribution of the neutral lattice. 

Setting Spf — Sfif = does not make 5(l)f and SFf 
zero outside the cluster. Nevertheless, we can write an 
expression for the polarization energy of ions that does 
not contain self-consistent potentials and fields outside 
the cluster. Subtracting Eq. (|l9| ) for the lattice with 
ions from the similar expression for the translationally- 
invariant lattice we obtain 



a i 

+ IT.T.^pT''^^'^''- (27) 



The first sum runs over all molecules a including the ions. 
The second sum over the ions appears because both p- 
and 0— components of the bilinear expression Eq. (|9|) are 

different at a = ion. The potentials 4>f'"^^ and fields i^'^'"' 
in the unrelaxed translationally-invariant lattice are eval- 



uated using Eqs. with p[ 



r,b - JO) 



and H° = 0. Thus, 



Eq. (g7|) gives the energy of a set of ions in an infinite 
lattice in which molecules beyond the imaginary cluster 
are not allowed to relax. The polarization energy of the 
set of ions in an infinite lattice is obtained as i? — > oo. 



A. P± in anthracene and PTCDA 



We start with a single ion, when Eq. ( [27[ ) yields either 
P+ or P_ in Eq. (|l|). The cluster of radius R is cen- 
tered on the unit cell that contains the anion or cation. 
Since clusters are defined in terms of unit cells, we know 
the number of molecules M{R) and Mvc/Nc = AttR^/S 
relates R to the molecular volume Vc/Nc in the crystal. 
The polarization energy P+ + P_ for ions at infinite sep- 
aration is shown in Fig. 3 as a function of M~^/^ for 
anthracene and PTCDA crystals. The "charges only" 
points refer to a = and polarization due entirely to 
charge redistribution; the other points are based on the 
B3LYP values of a (Table I). The largest clusters shown 
in Fig. 3 contain M — 2000 molecules, which corresponds 
to a cluster diameter of 2i? = 114 A for PTCDA and 97 
A for anthracene. 




Fig. 3 Convergence of P+ + P- for anthracene and PTCDA 
with M~^^^ , which is proportional to the inverse radius R^^ 
of the cluster. Straight lines are linear fits. Open symbols 
show the "charges only" results with a = 0. 

The polarization energy decreases with cluster size as 
more degrees of freedom for charge relaxation are added. 
At large R the missing part due to the molecules outside 
the cluster can be thought of as the polarization energy 
of a charge in the center of a cavity of radius i? in a con- 
tinous dielectric medium. Such energy is linear in 1/R. 



Linear extrapolation in Fig. 3 gives P+ -I- P_ 



-2.204 



eV for anthracene and -1.822 eV for PTCDA crystals. 
The smaller ion has the greater stabilization. 

The polarization energy of a charge in a spherical 
cavity in an anisotropic medium has been evaluated by 
Bounds and MunnEll 



P±(oo) -P±(i?) 



e 

2R 



1 - 



1 

Kcff 



(28) 



The effective dielectric constants is expressed in 
terms of the principal values ki < K2 < K3 of the di- 
electric tensor 



KcS ■ 



\/k2(k3 - Ki) 



F { arctan 



/v3 — Hi 



K3(k2-Ki) 
K2(rt3-Ki) 



(29) 



where F{(j),k) is the elliptic integral of the first kind. 
Equations (|28| ) and (^) determine the slope of the asym- 
totic behavior of P± in Fig. 3. As a consistency check we 
computed the slope using the dielectric tensor obtained 
in Sec. VIIL The values 2.482 eVA for anthracene and 
2.167 eVA for PTCDA are within 3% of the slope of the 
straight line in Fig. 3 drawn through the last two calcu- 
lated points. The "charges only" slopes are also within 
3% of the dielectric constants based on 3 = 0. 

Submolecules necessarily yield = P_, since ions 
are assumed to have equal and opposite charges and the 
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neutral lattice contains neither charges nor dipoles. The 
anthracene resultcil is P± = 1.42 eV for three points 
at the centers of rings and an effective a based on the 
static dielectric tensor of the crystal, while experimen- 
tal molecular a yields 1.51 eV for three submolecules. 
The electronic polarization, P+ + P_ ~ 2.8 eV, substan- 
tially exceeds our 2.20 eV. For PICDA, 11 submolecules 
and a similar to Table I lead t(M P+ + P- = 2.14 eV, 
which is again greater than our 1.82 eV. While the results 
clearly depend on the inputted a, the polarization energy 
of single ions found via charge redistribution is reduced 
compared to submolecules, especially when only one or 
a few are used. The difference is even greater when the 
are frozen at the gas-phase values. 
P- by 10 — 20 % in these systems. 



ion(Q) 



ion's charges p- 
This decreases P+ 
Although anthracene and PTCDA ions are at inversion 
centers, their atoms are not and thus experience local 
fields that are related by inversion. By contrast, sub- 
molecule charges are fixed at the outset and the field 
vanishes by symmetry at the center of the molecule. 

The transport gap, Eq. (|l|), of molecular crystals is di- 
rectly relatedo to photoelectron (PES) and inverse pho- 
toelectron spectra (IPES) on. surfaces, which yield adi- 
abatic P+ -I- P_ that includetl intramolecular relaxation, 
but not lattice relaxation. The inferredn3 P+ -|-P_ ~ 1.7 
eV for PTCDA films is quite consistent with the cal- 
culated 1.82 eV in the crystal. The importance of Et 
for electronic organic devices and recent thin film data 
were the motivation for the accurate calculation of elec- 
tronic polarization in the welLdefined limit of zero over- 
lap. The systems of interestEl have mobilities of 0.1 — 
1.0 cm^/Vs at room temperature, which is high for or- 
ganics and indicates that overlap corrections will have to 
be included. 

Partial charges and induced dipoles in the neutral lat- 
tice lead to P-|_ ^ P_. The individual components are 
shown in Fig. 4 and go as (X 1/R. The anion and 

cation slopes are equal at large R, in accord with Eq. (28). 
P_(_ and P_ are almost identical for anthracene and strik- 
ingly different for PTCDA. In the smallest cluster, which 
contains only the anion or cation, the ion interacts with 
the charges and dipoles Ji^ of the neutral lattice. 

Finite P+{M = 1) and P^{M = 1) are the energies 
of the cation and anion in the unrelaxed neutral lat- 
tice. They are nonzero due to fields in the neutral lat- 
tice. Without relaxation of the ion itself, P±(M — 1) 
is given by the second term of Eq. ( |27| ) with potentials 

= 4>i + 4>f^ . The relaxed ion in the field of the neu- 
tral lattice has charge distribution and P±{M = 1) 
is given by Eq. ( p7| ) with Sp"; — pf^* for a — ion and 
otherwise. We have assumed 5'°" = a for simplicity; 
more generally, nonzero A5 — a'°" — a will introduce a 
source term = AS^Fj in Eq. ( ^4b|) . 

The large PTCDA contributions at M = 1 do not can- 
cel exactly because the anion and cation charges are not 
precisely equal and opposite. Approximate electron-hole 
symmetry for the 7r-system of anthracene ensures almost 



equal and opposite charges. Our treatment gives two 
contributions to P± , an initial interaction at M = 1 that 
does not arise for submolecules and a relaxation or po- 
larization of the lattice that remains almost the same for 
the anion and cation. Such a distinction may be useful 
in future work. 

The charge distribution in PTCDA is such that posi- 
tive atomic charges are closer to the molecular centers. 
The resulting quadrupole moments of molecules create 
an average positive potential at each molecule in the neu- 
tral lattice. The 7r-electron density above and below the 
molecular plane also generates a quadrupole as discussed 
by Silinsh and Capek.l 



> 

CD 



ay 

I -1 



^ref^e^'"! , .,.8;; 



'^Vco^^'^ 



0.2 



0.4 



0.6 



0.8 



M 



■1/3 



Fig. 4 Convergence of P+ and P_ for anthracene and 
PTCDA with with M"^/^. Straight lines are linear fits. The 
values ai M — 1 are discussed in the text. 

In fact, the quadrupole contribution to P+ and P_ 
depends on the macroscopic shape of the sample. It is 
proportional to J d^r{l/r^), which gives finite contribu- 
tion from the remote parts of the sample. The contribu- 
tion to the potential is constant on the scale of the unit 
cell, because the corresponding contribution to the field, 
oc J (Pr{l/r'^), vanishes. Thus, the individual quantities 
P+ and P_ are not well defined, but the shape-dependent 
contribution cancels exactly in the sum, P_^ + P_, which 
enters Eq. (|l|) for the gap. 

The quantities / — P+ and A + P- can be viewed as the 
ionization potential and the electron affinity of the solid. 
We see that they depend on the macroscopic shape of 
the sample due to quadrupolar corrections. In general, 
the polarization energy of an abitrary set of charges in a 
crystal lattice depends on the shape of the macroscopic 
sample, unless the total charge is zero. 

The interpretation of P+ and P_ is, however, of con- 
siderable interest because PES and IPES spectra are 
relatedc3 to / — P+ and A + P_, respectively. Di- 
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rect comparisgn, therefore, requires a surface calculation. 
CompilationsEfEl of / — P+ for organics, while admittedly 
approximate, clearly point to > P_ and to the phys- 
ical meaning of individual polarization energies. 



B. Ion pairs in anthracene and PTCDA 

Polarization effects modify the interaction between 
charge carriers. We compute V{r) in Eq. (^) by replac- 
ing two neutral molecules in the lattice with a cation and 
anion. The cation is at the origin and the anion's cen- 
ter has crystallographic coordinates r = (a,b,c) given in 
Table II. As in the previous section, we consider an imag- 
inary cluster of radius R that contains both ions, solve 
the self-consistent Eqs. (23) and (^, and evaluate the 
energy E^pair using Eq. (27|). We repeat with larger R 
until V{r) converges. 
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Fig. 5 Interaction energy V^(r) 
pairs in clusters of M molecules. 



Eq. 



for various ion 



Figure 5 shows V{r) as a function of 1/M, which is 
proportional to the inverse cluster volume^ for various 
ion pairs. Since the pair is neutral, the 1/R contribution 
given by Eq. (28) vanishes, and the asymptotic behavior 
is linear in 1 / rF^. It represents the polarization energy of 
a dipole in the center of spherical cavity in a dielectric 
medium. 

The extrapolated values of V{r) for various pairs are 
presented in Table II, which also lists the distances be- 
tween centers and identifies pairs using the crystallo- 
graphic notation of Ref. The lowest charge-transfer 



(CT) exciton in PTCDA (Fig. 5) corresponds to neigh- 
bors in the stack. The next CT state is the second neigh- 
bor along the stack, which is closely followed in energy by 
other configurations for neighboring molecules in differ- 
ent stacks, as shown in Fig. 2. In anthracene the lowest 
CT state corresponds to-the closest neighbor. 

Bounds and SiebrandE3 obtained V{r) for anthracene 
using a single point per molecule and experimental po- 
larization data. Their 0.78 and 0.58 eV binding for the 
lowest CT states are less than our 0.922 and 0.821 eV; the 
large difference for the second neighbors is due to in-plane 
polarization and the inadequacy of a single point charge 
for long molecules that are end-to-end. Anisotropic a 
and charge redistribution in anthracene produce several 
instances (e.g. at 9.894 and 11.172 A ) where V{r) is not 
mono tonic in r. 

Table II Energies ^(r), in eV, of charge-transfer states 
within zero-overlap approximation for several separations r. 
The cation and anion are at (0, 0, 0) and (a, b, c) respectively. 



Pair (a, 6, c) 


r, A 


V{r) 


Kppr(r)'' 


Anthracene 








(i -i 0) 

(0,1,0) 


5.228 


-0.92 


-0.79 


6.016 


-0.82 


-0.76 


(1,0,0) 


8.553 


-0.55 


-0.50 


(1,0,1) 


9.458 


-0.51 


-0.50 


(- - 1) 

\ 2 ' 2 ' ^ 


9.894 


-0.54 


-0.57 


(-- ^ 0) 

\ 2 ' 2 ' 

(1,1,0) 


9.986 


-0.51 


-0.47 


10.456 


-0.49 


-0.42 


(0,0,1) 


11.172 


-0.51 


-0.56 


(1,-1,1) 


11.209 


-0.50 


-0.42 



PTCDA 

1,0,0) 
2,0,0) 



' 2 ' 



3,0,0) 
1 i -i) 

, 2 ' 2' 

0. 1.0, 

^' 2' 2^ 

1, -1,0) 

.-^) 

2, -1,0) 



2 i 

2 



3.726 


-1.06 


-0.75 


7.453 


-0.57 


-0.46 


10.558 


-0.46 


-0.48 


10.751 


-0.43 


-0.40 


11.179 


-0.38 


-0.32 


11.624 


-0.43 


-0.46 


11.998 


-0.42 


-0.46 


12.144 


-0.36 


-0.32 


12.563 


-0.39 


-0.42 


13.658 


-0.38 


-0.37 


14.124 


-0.33 


-0.32 



° Eq. (|30|) 

The point-charge approximatiorS gives almost 2.0 eV 
of binding for PTCDA neighbors in a-stack, twice the 
1.05 eV in Table II, while the binding is 0.99 eV for 
11 submolecules. Point charges are poor approximations 
for large molecules with interplanar separation of 3.4 A. 
Charge redistribution and partial charges provide a di- 
rect way to compute such electrostatic interactions. 

We may consider ^(r) for a cation-anion pair in a con- 
tinous anisotropic medium with dielectric tensor k. We 
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describe each ion with fixed partial charges pi and write 
an electrostatic expression in terms of the double sum 
over the atoms of the cation and anion. Since Ai^tot is 
defined relative to the neutral lattice, charge differences 
Ap-°" = — pj appear in V{t): 



Kppr(r) 



E 



(30) 



Here rf^ are the components of = — r~, and 
det{K) = K1K2K3 is the determinant of the dielectric 
tensor. Ej At large r, Eq. (|o|) reduces to point charges, 
since Jj^ — ensures vanishing interactions in the neu- 
tral lattice. The lowest-order corrections to Eq. ( ^0| ) are 
induced dipoles due to the other charge. 

Table II lists V'appr(r) values based on Eq. ( ^0|) and the 
dielectric tensors obtained in the next Section. Gas-phase 
charges and dielectric data are a simple and reasonably 
accurate approximation to the self-consistent calculation. 

We find comparable ^(r) ~ —0.5 cV for second neigh- 
bors in PTCDA stacks and first neighbors in different 
stacks (Fig. 2). Mazur and Petelenzll2l report instead 
that, by a small margin, an interstack neighbor has bind- 
ing 1.02 eV that even exceeds the first neighbor. They 
emphasize the competition between large in-plane polar- 
ization for neighbors in adjacent stacks and Coulomb in- 
teractions in the same stack. Such trade-offs are seen 
in Table II for both anthracene and PTCDA, although 
not as strongly as in Ref. |l^. They compute impair for 
fractional charges at 11 submolecules and then find the 
Coulomb interaction between the anion and cation using 
Lowdin charges in a 6-31G basis. While it is inconsis- 
tent to use different charges for polarization and direct 
interactions, it is natural to prefer atomic charges for 
the anion-cation interaction to arbitrarily placed partial 
charges. 

We also find (1,0,0) to be the lowest CT state by 
~ 0.3 eV on using, as in Ref. ^ 11 fractional charges 
and partial polarizabilities for both the polarization and 
cation-anion interaction. Hence (1,0,0) is the lowest CT 
state in PTCDA and, although important, the greater 
in-plane polarizability does not stabilize neighbors in dif- 
ferent chains below (2,0,0). 



zero overlap because space can be subdivided such that 
each molecule feels its own potential (/'(r). The charge 
distribution does not change when a constant is added 
to all (pi, according to Eq. (|l^). 

to the right-hand side of Eq. 



We add Ef; 



( |22b| ), and 

the corresponding term —rf"EQ to the right-hand side of 
Eq. ( ^2aD . We then solve ANNc self-consistent equations 
(^) and (^), and obtain the total dipole moment /i" of 
each molecule a in the unit cell using Eq. ( p!3[ ) . The total 
dipole moment of the unit cell fi = ^t" determines the 
polarization P = fJ,/vc, where Vc is the unit-cell volume. 
Repeating this procedure with Eq directed along each of 
the coordinate axes we find the tensor ( that relates P 
to Eq: 



P = C^^o. 



(31) 



Alternatively, one may choose to differentiate the self- 
consistent equations with respect to Eq to find an ex- 
plicit expression for C. 

The susceptibility tensor x is defined by the relation 
P = xF, where F is the total average macroscopic field 
that is creaied by external sources as well as the polarized 
solid itself.tJ According to the Lorentz relation, a dielec- 
tric sphere with uniform polarization P creates a field 
— 47rP/3 in its center. The spherical shape is consistent 
with the definition of V(r) in Eq. (|l|). Thus, 



An 



(32) 



Using Eq. (|3^) we eliminate Eq in favor of F and using 
K = 1 + Attx wc obtain the dielectric tensor 



1 + (8^/3)C 
1 - (47r/3)C 



(33) 



Table III Experimental (averaged) and calculated compo- 
nents of dielectric tensor of anthracene and PTCDA 



Inputs 



Kl 



K3 



AnthrE 

ExperimcntBtrna 2.49(10) 3.07(10) 4.04(20) 28(2) 
B3LYP/6-311++G** 2.23 2.91 4.03 31.6 
Charges only (5 ^ 0) 1.36 2.39 3.90 34.5 



VII. DIELECTRIC TENSOR 

We now summarize the calculation of the dielectric ten- 
sor by considering a sample in a uniform electric field. 
These results have been reported previously.tll The pro- 
cedure is similar to the polarization energy of the neutral 
lattice in Sec. V. 

Formally, an applied field breaks translational invari- 
ance, since the electrostatic potential is unbound. Never- 
theless, we can add an appropriate constant to the poten- 
tial in each unit cell and restore translational symmetry, 
without effect on charge distribution. This follows for 



Experimentlij 

B3LYP/6-311-K-fG** 
Charges only (S — 0) 



1.9(1) 
1.85 
1.96 
1.01 



4.3(2) 
4.07 
3.98 
3.74 



4.6(2) 
4.07 
4.02 
3.81 



Table III lists the principal components of k in an- 
thracene crystals and PTCDA films. K2 is along the b 
axis; 9 is the angle m.±iie_a,c plane between ki and a. 
The anthracene datalljcJta are averages over indepen- 
dent measurements of the dielectric tensor and refrac- 
tive indices. The calculated values are for a in Table I, 
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and additional a inputs for larger bases are reported in 
Ref. |ll|. The results agree with the available experimen- 
tal data, which has an accuracy of a few percent. While 
3 = accounts for K3, the largest component, the quan- 
titative importance of accurate a is clearly seen in Ta- 
bic III. 

Equation (^3|) for k is strictly based on neutral 
molecules in a translationally-invariant crystal. The en- 
ergies of molecular ions were found instead within clus- 
ters of radius R. As noted above, Eq. ( |2^ ) reproduces the 
giope in Figs. 3 and 4 within 3%. This demon- 
strates the internal consistency of the procedure. 



VIII. DISCUSSION 

In the limit of zero overlap, molecules in organic 
crystals are quantum systems with purely classical 
electrostatic interactions. We have developed a self- 
consistent approach that treats each molecule quantum- 
mechanically, subject to the external fields of all other 
molecules. We have found that such external fields can 
be treated perturbatively to a good accuracy, expanding 
the solution of Schrodinger equation for each molecule 
near the gas-phase solution. Self-consistent analysis of 
large systems, with over 10^ atoms, is straightforward. 

The self-consistent procedure captures an important 
effect: the redistribution of charge in molecules subject 
to (nonuniform) external fields. Direct description of 
charge redistribution avoids the ambiguity that accompa- 
nies partitioning the molecular polarizability over a num- 
ber of polarizable points. 

As outlined in Section III, we discretize the molecular 
charge distribution p(r) by assigning a partial charge pi 
and induced dipole fii to every atom i. The atom-atom 
polarizability tensor 11^^ then governs charge redistribu- 
tion in external fields. The tensor 11^ is conveniently 
found quantum-mechanically at a semiempirical level, 
such as INDO/S, which is well suited for introducing po- 
tentials or site energies. The correction a = a — a'^' 
distributed over the atoms, and as in previous theory, 
we still face the familiar problems associated with sub- 
molecules. The present approach, however, allows for a 
certain compensation of error: the partitioning now only 
involves the small correction a, which is about 10 — 20% 
of the actual molecular polarizability a in large conju- 
gated molecules. 

The electrostatic potential (/)(r) created by a molecule 
with charge distribution p(r) is given by Eq. (|). While 
atomic charges cannot be defined uniquely, charges or 
dipoles that reproduce (f)(r) accurately at intermolecular 
distances in crystals would be completely satisfactory. 
Thus, ab-initio p{r) geperate (/)(r) outside the molecule, 
which provides a basiaHj for assigning discrete charges. 
This suggests an interesting possibility of introducing 
gas-phase atomic dipoles fif ^ along with the gas-phase 
atomic charges pf'^ in Eq. (8b), as long as they improve 



the description of fields created by a neutral molecule. 

Another possible extention is to introduce atomic 
quadrupoles to represent 7r-electron density above and 
below the plane. This will result in an additional lOiV 
scalar equations per molecule. In practice, only the qnn 
component normal to the conjugation plane is needed for 
TT-electrons. Classical multipoles lead to complications 
even in the limit of no overlap. These are corrections, 
however, to charges and induced dipoles whose fields Fi 
we have found, and perturbation theory may well be suf- 
ficient. 

We have already contrasted our method to the existing 
approach in which a molecule is represented by a number 
of polarizable points, termed submolecules. The qualita- 
tive differences are the electrostatic energy of the neutral 
lattice and the energy P±{M = 1) of an ion in the un- 
relaxed lattice in Fig. 4. Charge redistribution avoids 
entirely the number and location of the submolecules. 
Since the actual a enters in either case, the numerical 
results for P+ + P_ in Section VI and for ^(r) in Ta- 
ble II do not differ greatly, especially in comparison to 
calculations with many submolecules. 

Crystalline organic films that function as electron or 
hole conductors are of particular interest for organic 
electronics Jji3 Although overlap is then finite, it should 
be considered as a correction to much greater polarization 
energies. We have a systematic method for computing 
the electronic part of P+-f P_ , while the total polarization 
energy in Eq. (^ is about 10% larger due to the lattice 
contributions, even before overlap or charge transport is 
introduced. Reliable comparisons of the large electronic 
part are the first step. We comment on two issues. 

First, the greater stabilization of separated ions has 
major implications on the binding energy |y(r)| of neigh- 
boring kms. Since the energy of the lowest CT in 
PTCDAlj is comparable to the lowest Frenkel exci- 
ton, \V[r)\ is closely related to discussions and debates 
about the magnitude of exciton binding energies .Eil Lat- 
tice (molecular) relaxation about a separated anion and 
cation is probably in the range of 200 meV each, while 
relaxation about an adjacent ion pair is less than 100 
meV; the lattice contribution decreases |T^(r)|. Finite 
overlaps and transfer integrals in the 50 — 100 meV range 
are typical for hopping transport and stabilize polarons 
in Holstein models. CT states, by contrast, have vanish- 
ing bandwidth or mobility because two-electron transfers 
are needed. Overlap also preferentially stabilizes individ- 
ual ions and hence reduces |V^(r)|. The lowest CT state 
of PTCDA or anthracene has electronic binding of '--^ 1 
eV in Table II that could be reduced substantially due 
to lattice relaxation or by overlap in systems with mobile 
charges. 

The second issue addresses novel features of self- 
consistent atomic charges. As noted above, charge re- 
distribution of isolated ions increases P_|_ or P_ by about 
10% even when ions are at inversion centers. The lower 
symmetry of ion pairs leads to more extensive charge 
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redistribution of mutual polarization that is automat- 
ically included in our treatment through Ily for ions. 
Charge redistribution on ions is important. Its contribu- 
tion can be estimated by self-consistent calculation with 
ionic charges and dipoles set to gas-phase values. Specific 
contributions are found by comparing two self-consistent 
solutions. 

In summary, we have implemented a general approach 
to electronic polarization of organic molecular crystals 
in the limit of zero intermolecular overlap. Redistribu- 
tion of partial atomic charges is governed by the atom- 
atom polarizability tensor Iljj , which provides a quantum 
mechanical basis for the description of electronic polar- 
ization. Partial charges replace the microelectronics of 
postulated submolecules in previous treatments and in- 
troduce such new features as electronic stabilization of 
the neutral lattice and charge relaxation on ions. Self- 
consistent atomic charges and induced dipoles relate the 
molecular polarizability of anthracene or PTCDA to the 
dielectric tensor and the encrgic;s of fixed ions and ion 
pairs in the crystal. The cluster approach used for ions 
is suitable for surfaces and other systems with reduced 
symmetry, provided that all molecular positions are spec- 
ified. Zero overlap provides a starting point for the treat- 
ment of electronic polarization in organic systems with 
mobile localized charge carriers. 
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